Relationships between nonmesonic-weak-decays in different hypernuclei 



Franjo Krmpotic 
Institute de Fisica La Plata, CONICET, 1900 La Plata, 
Argentina, and Facultad de Ciencias Astronomicas y Geofisicas, 
Universidad Nacional de La Plata, 1900 La Plata, Argentina. 
(Dated: November 18, 2010) 

Using as a tool the s-wave approximation (sWA), this work demonstrates that the nonmesonic 
weak decay transition rates T n and T p can be expressed in all hypernuclei up to A 9 Si (and very likely 
in heavier ones too) in the same way as in the s-shell hypernuclei, i.e., as a linear combination of only 
three elementary transition rates. This finding leads to the analytic prediction that, independently 
of the transition mechanism, all hypernuclei that are on the stability line (N = Z), i.e., ^He, ^Li, 
\Be, j^B, A 3 C, a 7 0, i 9 Si, e ^ c should roughly have the same ratio T n /T p , the magnitude of which 
rapidly increases when one approaches the neutron drip-line (N ^> Z), and opposite happens when 
one goes toward the proton drip- line (N < Z). 



I. INTRODUCTION 

The nonmesonic weak decay (NMWD) of A hypernu- 
clei, AN — > nN, takes place only within a nuclear envi- 
ronment with the decay rate Tn (N = p,n). Without 
producing any additional on-shell particle (as does the 
mesonic weak decay A — > irN), the mass is changed by 
176 MeV, and the strangeness by | AS| = 1, which implies 
the most radical modification of an elementary particle 
within the nucleus. At the same time it offers the best 
opportunity to study the strangeness-changing interac- 
tion between hadrons, and is the main decay channel for 
medium and heavy hypernuclei. 

With the incorporation of strangeness, the radioac- 
tivity domain is extended to three dimensions (IM,Z,S), 
which, because of the additional binding due to the A- 
hyperon, is even richer in elements than the ordinary 
(N, Z) domain. (For instance, while the one-neutron sep- 
aration energy in 20 C is 1.01 MeV, it is 1.63 MeV in 
[l|.) This attribute of hypernuclei has motivated a recent 
proposal to produce neutron rich A-hypernuclei at the 
Japan Proton Accelerator Research Complex (J-PARC), 
including A He Q. 

Important experimental efforts have been invested in 
hypernuclear weak physics during the last few years 
Il6| . The correlative theoretical advances in our knowl- 

te of the NMWD, have been also quite significant [l7r - 
For recent review articles see Refs. |50H53j . The 
ratio r„/ p = r„/r p , together with the asymmetry pa- 
rameter ratio <2a for emission of protons from polarized 
hypernuclei [13, HH, has been in the past and still are the 
main concerns in the physics of NMWD. For a long time 
the large experimental value for the T n / p ratio (close to 
unity) remained unexplained. But, recent improved data 
tend to converge to = 0.5 0-G3 , both for A He (s-shell) 
and A 2 C (p-shell), indicating similarity in the transition 
mechanism. 



It is also speculated that the NMWD could have an important 
role in the stability of rotating neutron stars with respect to 
gravitational wave emission |j| . 



In the meantime the theoretical estimates of T n j p , done 
within the one meson-xchange (OME) model, have in- 
creased. For instance, Parrefio, and Ramos [24j have 
found r„ /p ( A He) = 0.34 - 0.46, and r„ /p ( A 2 C) = 0.29 - 
0.34, when the exchanges of the complete pseudoscalar 
(•7r, K, rj) and vector (p, lu,K*) meson octets are taken 
into account, with the weak coupling constants obtained 
from soft meson theorems and SU(Q)w [13; EH- The 
dominant role is played by the exchange of pion and 
kaon mesons, and when their effect is combined with 
the Direct-Quark (DQ) model, the value of the n/p ratio 
is increased up to 0.70 [20, [23|. However, these transi- 
tion mechanisms continue to predict too large and neg- 
ative value for a\. There are two recent proposals to 
bring this value into agreement with experiments by go- 
ing beyond the OME+DQ models. The first considers 
incorporating new scalar-isoscalar terms induced by 27T- 
exchanges [38[.( See also Refs. [H,|3l| on the relevance of 
these terms.) In the second, in addition to the model of 
7r + 2tt I p + 2n I 'a + 2p/ 's + uj + K exc hang es, is introduced 
the axial- vector ai-meson exchange [40(. 

Quite recently we have discussed the parameter cla 
within the independent particle shell model (IPSM), to- 
gether with the s-wave approximation (sWA) [37]. The 
corollary of this study was that, independently of the 
NMWD dynamics, this observable has the same value 
in all hypernuclei that have totally full proton subshells, 
such as A He and A 2 C, and very likely also in the remain- 
ing hypernuclei. This result is a direct consequence of the 
fact that a \, same as T n / p , is a ratio of two transition 
rates, which makes it, in absence of Final State Interac- 
tions (FSI), dependent purely on the dynamical features 
of the NMWD. 

The aim of this work is twofold. First, we establish the 
link between the theoretical formalism for the NMWD 
of the s-shell hypernuclei originally introduced by Block 
and Dalitz [55J , and the general formalism used presently 
for any type of hypernuclei. Second, we show that the 
IPSM framework, together with the sWA, allow us to for- 
mulate the rates T n within the p, d, etc shells in terms 
of the s-shell nuclear matrix elements (NME). Previous 
research in this direction has been done by Alberico, and 
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Garbarino [2l| and by Cohen |56| . Later on, it is demon- 
strated that regardless of the decay mechanism: i) all 
hypernuclei with the same number of protons and neu- 
trons (i.e., with Z = N) should have the same ratio T n i p , 
ii) the value of this observable increases (decreases) as 
the neutron (proton) excess is enlarged, and iii) simple 
analytic relationships exist between r„, r p , and Y n / p in 
different hypernuclei with the same mass number A. The 
derivation of these results, same as those on the parame- 
ter a a 37J , is based on the assumption that the emission 
of the nucleons N from different single-particle states is 
affected in a similar way by the FSI. Then, before pre- 
senting the formalism, it might be convenient to comment 
on the relationship between the ratio T n /T p and the FSI. 

The primary partial decay rates Tjv are in principle 
derivable from the measurements of emitted nucleons n 
and N spectra. These are: i) the single-nucleon spec- 
tra S{Eff), as a function of one-nucleon kinetic ener- 
gies En, and ii) nN coincidence spectra S(E n tf), and 
S (cos n N), as functions of the sum of kinetic energies 
E n N = E n + En, and the opening angle 9 n N, respec- 
tively. From these spectra are determined the numbers 
of protons N p , and neutrons N„, and numbers of pairs 
N nn , and N np , which are not related in a simple way 
with r„ and T p . This is because not all primary nucle- 
ons, originated by the NMWD, are measured. In propa- 
gating within the nuclear environment they interact with 
the surrounding nucleons, and in some cases they change 
their momenta, and energies, some of them even can 
be absorbed by the medium, and emission of additional 
(secondary) nucleons can take place as well [43Tj48j . All 
these processes represent a complicated many-body prob- 
lem, and are generically designated as FSI. To describe 
them, keeping the calculations feasible, are indispensable 
model assumptions, and the FSI are usually simulated by 
a semi-classical model, developed by Ramos et al. [19j, 
and called the intranuclear cascade (INC) model. This 
model interrelates the rates T n , and T p with the numbers 
N n , N p , N, m , and N np , and therefore, as stressed recently 
by Bauer and Garbarino ^3], the FSI described by the 
INC model should not be included in the evaluation of 
decay rates T n , and T p . However, not all FSI are con- 
sidered within the INC model, and which additional FSI 
contribute to the NMWD spectra and decay rates, and 
how and which of them should be included in the calcu- 
lation are nontrivial questions. Some candidates are as 
follows: 

(i) Short range correlations (SRCs) acting on final nN 
states; here one starts from the plane wave approximation 
for the outgoing nucleons and the SRCs are incorporated 
a posteriori, either phenomenologically through Jastrow- 
like SRC functions, or by solving the Bethe-Goldstonc 
equation. The first approach is used within both nuclear 
matter [3(1 I43T - |45| . and finite nuclei calculations [iH 
[1^-[2^|, and the second one only in the shell-model- type 
calculations [H, [H El S3] ■ 

(ii) Self-energy and vertex particle-hole corrections, 
and RPA-like rescattering effects (see, for instance, the 



diagrams (b)-(d) in Fig. 2 of Ref. [5J|). It is not known 
whether these FSI contribute coherently or incoherently, 
and it can even happen that (b) and (c) cancel out, as 
do the divergences in the vertex, and fermion self-energy 
corrections in the QED, because of the Ward identity. 
(Something similar happens also in the nuclear particle- 
phonon-coupling model.) The first ones can be associated 
with the mean-field effects on the single-particle wave 
functions engendered by an energy-dependent complex 
optical potential [49| . 

(iii) Interactions of the deep-hole states (which become 
highly excited states in the continuum after the NMWD) 
with more complicated configurations (2h\p, 3h2p ■ ■ ■, 
collective states, etc), which spread their transition 
strengths in relatively large energy intervals (4lj . 

There is no theoretical study in the literature on the 
NMWD encompassing all aspects of the FSI. The devel- 
opment of a microscopic many-body model for the FSI 
described by the INC model would be also extremely wel- 
come, and so far only in Refs. [H,|4^] were the first steps 
taken toward this goal. Finally, the two-body induced 
NMWD ANN -> nNN, which has been recently mea- 
sured [HI , should be also considered 0. Briefly, the is- 
sue of FSI is a tough nut to crack, and a lot of theo- 
retical work has to be done still, particularly in relation 
to the recentl y m easured spectra S(En), S(E u n), and 
S(cos9 n N) [a-tla]. which are certainly affected by them. 
However, as the purpose of the present contribution is 
not to make progress in this direction, among all possi- 
ble FSI, only the SRCs will be considered here. This is 
done phenomenologically, and initial A7V state SRCs are 
included on the same footing [3, [26l - [28| . It is our belief 
that this is a fair approximation for the objectives of the 
present work. 



II. DECAY RATES 

To derive the NMWD rate within the IPSM we start 
with the Fermi Golden Rule. For a hypernucleus, in its 
ground state with spin Jj and energy Ej 1 , decaying to: i) 
several states a at in the residual nuclei with spins Jf and 
energies E aN j F , and ii) two free nucleons n and N, with 
momenta p n , and pn, kinetic energies E n = /2M, and 
En = p%/2M, and total spin S, reads |26l428j 

T N = 2tt f S ( A aN j F - E R - E n - E N ) 

O 7_ <J 



Soln Jf 

1/ G T IT/I T \|2 ^P" 

x \(PnPNS;a N J F \V\Ji)\ 



(2tt) 3 (2tt) 3 ' 



(2.1) 



where for sake of simplicity we have suppressed the mag- 
netic quantum numbers. The NMWD dynamics, con- 
tained within the weak hypernuclear transition potential 



2 It also has been shown that the kinematical and nonlocal and 
kinematical effects on the NMWD could be sizable [27]]. 
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V, will be described by the OME model. The wave func- 
tions for the kets \p n p N SM s J F M F ) and | J/M/) are as- 
sumed to be antisymmetrized and normalized, and the 
two emitted nucleons n and N are described by plane 
waves. Initial and final SRCs are included phenomeno- 
logically at a Jastrow-like level, while the finite nucleon 
size effects at the interaction vertices are gauged by 
monopole form factors [18j, l26j. Moreover, 



Er 



\Pr, 



Pn\ 



E n + En + 2 cos OunV E u En 



2M{A-2) 



A 



is the recoil energy of the residual nucleus, and 

A aN J F = A + Ej T - E aN j F , 



(2.2) 



(2.3) 



with A = Ma — M = 176 MeV is the liberated energy. 

It could be convenient to perform a transformation to 
the relative and cm. i) momenta: p = (p n — Piv)/2, P = 
Pn +Pn, ii) coordinates r = r n - r N , R = (r„ + r N )/2, 
and iii) orbital angular momenta I and L. The energy 
conservation is expressed as 



E„ 



E, 



where 



El 

M 
P 2 A 
4M~ 



E R = 



2M(A - 2) ' 
A 



-Er 



(2.5) 



are, respectively, the energies of the relative motion of the 
outgoing pair, of the recoil, and of the total cm. motion 
(including the recoil). 

Following the analytical developments done in 
Ref. [26j], the transition rate can be expressed as a func- 
tion of the cm. momentum P: 



a N J F 



X P\ P 



P 2 T aN j F {pP), (2, 



A QjvJf = e p + e P - A aN j F = 0, (2.4) with 



SMLTJ 



Y J M( P lPLXSJT;j A j N Jt AN )(J I \\ (oJ w o} A ) \\a N J F ) 



JN 



(2.7) 



where 



and 



= 2 



A-2 



CtN Jf 1 



A 



A-2 



P 2 



-P 2 



(2.8) 

(2.9) 
has to be 



It is clear that the condition P < Pa q 
fulfilled for each final state |ajv Jf)- Moreover 

M( P lPLXSJT;j AjN Jt AN ) = -±= [1- (-) l + s + T ] 



x L (P)(lL\SJT\V(p)\j AjN Jt AN ), 



(2.10) 



where (and henceforth) the ket |), unlike |), indicates that 
the state is not antisymmetrized, 



L (P) = J R 2 dRj L {PR)R 0L (b/V2,R), 



(2.11) 



is the overlap of the cm. radial wave functions Rol, and 
ji, for the bound and outgoing particles, respectively, and 
b is the harmonic oscillator size parameter. More, A = 
I + L, T = {TM T ,M T = m tA + m tN }, and t A jv = {*a = 



1/2, m tA = -1/2,% = 1/2, m tN }, with m tp = 1/2, and 
m tn = —1/2, where we have assumed that the AN — > nN 
interaction occurs with the isospin change AT = 1/2. 
Explicitly, 



|tA 



\T = 1), for N = n 

(|T = l)-|T = 0))/v / 2, for N=p ' 

(2.12) 

It might be pertinent to mention that the factor (A— 2)/ A 
in Eqs. (|2.6p . (|2.8|) . and p.9|) comes from the recoil effect, 
which, in the same way as the spreading of the deep hole 
states, is relevant for the NMWD spectra [IHH^. but its 
role is of minor importance for the total transition rates 
Tat. 



A. Independent-article shell model 

Up to now nothing has been said about the initial state 
| J/), and final states \cxnJf)- Within the IPSM the fol- 
lowing assumptions are made, which greatly simplify the 
numerical calculations: 

1) The initial hypernuclear state is taken as a A- particle 
in single- particle state j A = 0s 1 / 2 weakly coupled to an 
(A — 1) nuclear core of spin Jc, i.e., \ Jj) = \{Jcjh)Ji)- 
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2) When the nucleon inducing the decay is the single- 
particle state j'jv (j = nlj), the final residual nucleus 
states are: IoatJf) = ^JcJm^Jf)- 

3. We adopt the simplest version of the IPSM, in which 
all the relevant particle states are assumed to be station- 
ary, and the liberated energy is 



A 3 

N 



-JV> 



(2.13) 



where N = p,n, and the e's are single-particle energies. 
(The non-stationary version of the IPSM is discussed in 
Ref. S3.) 

Within this scheme, we get [HI, [28[ 



■ N 



E^ 



1 N 



J=j+l/2 

E n^ij, (2.i4) 

J=|j-1/2| 



where the summation goes over all single-particle transi- 
tion rates T J Nl which in turn results from the sum over 
the values of J = jn + j A of products of the spectro- 
scopic factors FLj with the partial An nN transition 
rates Ti- 3 NJ - For the s-shell nuclei the later have the same 
physical meaning as the quantities Rnj introduced in the 
seminal work of Block and Dalitz [55| (see also Ref. [5g| ) , 
i.e., K%'f = R NJ . E 

The spectroscopic factors F NJ are defined as 



*ij = ^ E K^H («) r W J f)\ 2 (2-15) 



with the notation J = V2J+T, while the partial transi- 
tion rates read 



TZ 3 - 



2M,s 



.4 



dPP 2 J(P 3 f ) 2 - P 2 



A-2, 

E \M(plPL\SJT;j AjN Jt AN )\ 2 



SILXT 



with 



P J — 9 

N 



A-2 
A 



M N A% 



the maximum value of P for each j at , and 



V 



A 



A-2 



(P 3 N ) 2 



p2 



(2.16) 



(2.17) 



(2.18) 



the corresponding relative momentum. 

It should be stressed that the most important virtue 
of the IPSM is that the index a becomes superfluous, 



3 In order to use here the same notation for Rnj as in Ref. |55|l . as 
well as to write TZ J N j instead of Ti J j N , the jjv variable employed 
in previous publications is frequently split here in j, and N. 



and the summation on the final spins Jp can be carried 
out without knowing the nuclear structure of the initial 
and final nuclear states. This simplifies enormously the 
numerical calculations. As far as we know, the IPSM has 
been used to a great extent in all previous finite nucleus 
evaluation of the NMWD. 



B. s-wave approximation 

Galeao [57[ has shown that the matrix elements in 
(|2.10j) can be cast in the form 

(lLXSJT\V(p)\j A j N Jt AN ) = J2( lSKT W(p)\\S'Kt AN ) 

KS'\ 



x Q(lLXSJl N j N ;KS'), 



(2.19) 



with 



Q(ILXSJI NJN ;KS') = 

J 3N I 



S In 



x i N \j N s'k 2 ((j\oLi N \oooi N i N ) I i ^ - 

LI X \ ( L K J 
S J K MS' l N I 



(2.20) 



where (0 • • • | • • • In) are the Moshinsky brackets [So ]. The 
plain sWA implies that we make I = in (|2.19[) . and 
jOni), which leads to 

(lLXSJT\V(p)\j A j N Jt AN )= E C (lXSJl N j N ;K), 

if =0,1 

x (lSKT\V( P )\0KKt AN ), 

(2.21) 

with 

C (lXSJl N j N ;K) = K 2 Xj N {000l N lN\000l N lN) 





3N 


* ) 




I 




it 


K 


In J 




J 


In J 



(2.22) 

In particular, for the s-shell hypernuclei 

CqQISJO, 1/2; K) = 5jk, (2.23) 
and using Eqs. ([2~10| . and ([2T2]) . Eq. (|2~16]) becomes 



'^NJ = K NJ = [1 + O n N) Y 



A-2 



X 



dPP 2 ^P 2 -Pi@ 2 (P) (2.24) 
^[1- (-y+^]\(lSJT\V(p)\0JJT)\ 2 , 



SIT 
'N = r s 



where Pm = Pf r 1/2 , and 



3 "( P )-(f) 1/2 ^ (Pfc)2/2 



(2.25) 
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As it is well known, the corresponding transition rates 





3 

^-Rivo - 


f- -i?JVl, 




r n (iH) = 


2^n0 " 


- 2^ n1 ' 


Tp(aH) — Rpo, 


r p (lHe) = 


2^p0 "+ 


■ \RpU 


T„( A He) = Rno, 


I>(iHe) = 


- 


\- -^Rni, 


(2.26) 



with N = n,p, depend only on four single-particle tran- 
sition rates R n o, R n i, Rpo and R p \. 

Here we will express the Fat's of heavier hypernuclei in 
the same way as was done in Eq. (|2.26|) for the s-shell 
hypernuclei, i.e., as a linear combination of Rno and Rni 
only: 



Tn — FnoRno + J~niR 



Nl- 



(2.27) 



To derive the generalized spectroscopic factors (GSF's) 
Tno and Fni for hypernuclei up to A 9 Si we perform sum- 
mations over A in (|2.16[) for each single-particle state 
3 = P3/2,Pi/2, and j = d 5/2 . The resulting Tl° NJ turn 
out to be quite similar to (j2.24|) for TZ^'J , except that 
now Pjv, and 0o(P) are substituted, respectively, by P 3 N , 
and 0i N (P). Thus, we supplement the plain sWA with 
the substitution, 



N ■ 



(2.28) 



which are fair approximations for the evaluations of ratios 
r„/r p and oa- □ In this way we get 



K 



P3/2 


Rno 


Rni 


Nl ~ 


3 


6 ' 


Pl/2 


-Rao 


Rni 


Nl — 


6 


3 ' 


d 5/2 


3Rno 


Rni 


N2 ~ 


20 


10 



T?P3/2 _ 

'*-A2 — 



,X -N0 ~ 



/V W3 



R 



Nl 



2 ' 

Rni 
2 ' 
Rni 



(2.29) 



where the numerical factors come from the summation on 
A of the squares of coefficients Cq(IXSJInJn', K) given by 
(12321 . 

It could be useful to exp ress the Eq. (|2.24l) within the 
Block-Dalitz notation |55| : 

a = ( 1 So|^| 1 So), b = ( 3 P |F| 1 So>, c=(^ 1 |F|% 1 ), 

d = ( 3 D 1 |^| 3 S 1 ), e-( 1 P 1 |F| 3 S 1 ), f = ( 3 P 1 |y|3S 1 ), 

(2.30) 



4 This sWA has been used in Ref. [37| to relate the matrix ele- 
ments M(plPLXSJT;jAjN = 0p 3 / 2 , JtAjv), and (plSJT\V\\ = 
0, JJtAjv) in A 2 C (see [37|, Eq. (B2)]. There are two misprints 
in [37|, Eq. (B2)]. The correct results are: M(p2, PI, 1110; Ap) = 
J i f d(p)(Pl|ll), and M(p2,Pl,2120;Ap) = ^d(p)(Pl|ll). 



for the NME. Assuming the same value of single-particle 
energies for protons and neutrons in the state 0siy 2 , one 
gets the well known results 



Rpo = a, 2 + b 2 , 



R n0 = 2(a 2 + b 2 ), " 2 

Rni = 2/ 2 , Rpx =c 2 + d 2 + e 2 + f, (2.31) 



where 



A 



A-2 



7T 

v / 

dPP 2 ^JP 2 -P^0 2 (P) a 2 (p), (2.32) 

and similarly for &,•••/. As the NME depend very 
weakly on the momentum p we can compute them at 
p = pa = \/MnA [13] , and write 

a = J a 2 (pA), etc 

with 

2M N 



(2.33) 



-Pn 



A 
A-2 



N dPP 2 ^/p 2 -Pi0l(P), 
which after performing the integration reads [5£ 



J,, = -2M( \ 2ir^-j^A N be- z h{z), 



(2.34) 



(2.35) 



where Ii(z) is modified Bessel function of the first kind, 
and 



z = A N M N b 2 



A-2 
A 



(2.36) 



Using the asymptotic form of Ii(z) (Ii(z) = e z /^/2ttz) 
one gets 



Jo 



2M N p A , 



(2.37) 



which is the same result as that derived previously in 
Ref. 37[, where the recoil effect was not considered. 



Therefore one sees that after performing the integration 
in (|2.32l) this effect is totally washed out from the tran- 
sition rates. 

To evaluate the IV from (|2.14p . as well as to derive 
the GSF's in (|2.27[) . we need to know the spectroscopic 
factors F N j given by (|2.15[) . These, in turn, depend on 
the angular momenta Jc and Jj, which are fixed from 
the experimental data, and are exhibited in Table HI The 
F NJ - values for most of the hypernuclei discussed here 
are listed in (2|| Table I]. The remaining can be easily 
inferred from this table, except for the p-wave ones in A Li 
(F P N T = 5/8, and F P N % 2 = 3/8), and in A Be (F N a ( 2 = 
3/4, and F^T = 5/4). The resulting GSF's F NJ are 
listed in Table Q] It is noticeable that in the jj closed 
shells the singlet to triplet ratio (J^no '■ -Pzvi) is always 
(1:3). 
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TABLE I. Core spins Jc, initial spins Ji, and the generalized 
spectroscopic factors J-nj within the jj-coupling. 



A ^ 


In It 


*/ n(J 


7T . 

•> nl 


T n 


T , 
•> pi 


5 «He 

A™ 


n i 

u, 2 


1 /2 

1/ z 


3/2 


1 11 

1/ z 


3/2 


I Li 


1 I 

!' 2 


17/24 


43/24 


17/24 


43/24 




i 

u, 2 


1 


3 


1 12 

1/ z 


3/2 


?Re 


n i 

u, 2 


5 /8 


15 /8 


5 /8 


15/8 




1 

U, 2 


1 /2 


3/2 


1 


3 


11 R 
A D 


3 5 

O, 2 


25 /24 


TO /94 


25 /24 


TO /94 


A ^ 


^ 1 

2 ' ± 


13/12 


29/12 


1 


3 


A ^ 


0,| 


1 


3 


1 


3 


21 r< 
A ^ 


0,5 


13/8 


39/8 


1 


3 


A 6 




7/6 


10/3 


5/4 


15/4 


A 7 


0,| 


5/4 


15/4 


5/4 


15/4 


fSi 


2> Z 


33/20 


23/5 


13/8 


39/8 


fSi 


0,| 


13/8 


39/8 


13/8 


39/8 



These and similar results for other values of A are very 
likely independent of both the decay mechanism and the 
final state interactions. 

Detailed calculations of Refs. [IH HE] have proved 
that the contribution of the p partial wave to NMWD in 
p-shell hypernuclei, as well as in heavy-mass systems, is 
relatively small (^10%). In particular, Itonaga et al. (2f| 
have explored the decay rates T ni and T p in hypernuclei 
from A = 4 up to A = 209, establishing that the p- 
wave contributions to the calculated total one-nucleon 
induced decay rates r nm = T p + T n , and ratios T n / p are 
only a few percent of the respective s-wave contributions 
(see [H, Fig. 9]). They have attributed this finding to 
the short range of the decay interaction. In fact, it is well 
known that the ranges of the radial pieces v(r) of the 
OME potentials are inversely proportional to the meson 
masses (see [25|, Figs. 3, 4, and 5]), being the largest 
inverse mass that of the pion (m^ 1 = 1.4 fm). Then, 
when one analyzes the radial matrix element in (|2.19p . 
which reads (see [13, Eq. (A19)]) 



In all hypernuclei heavier than ^Hc the contribution 



of the state 0si/ 2 to the total rate Tjy is 



2 Rm ' 



(2.38) 



while the contributions of the single-particle states 0p 3 / 2 , 
0pi/2i etc depend on their occupations, which, in turn, 
are reflected in the values of the GSF's listed in Table |U 



For instance, T p 3/2 



r;, and Tn /2 = in all hyper- 
nuclei with Z > 6, and N > 6, respectively. In the same 



way, F 



Pl/2 
P 



r;/2, 



and r 



Pl/2 



in all hypernuclei 



with Z > 8, and N > 8, respectively. The orbital 0d 5 / 2 
supplies less transition strength than the 0si/2 state, and 
it is given by 



1 N 



- (Rno 



3i? 



Nil 



(2.39) 



in all hypernuclei with Z > 14, or N > 14. 

Several years ago, by means of the sWA, Cohen [Hf 
arrived at the estimate 



Fn/p 



RnO + 3i?„i 



(2.40) 



Rpo + 3i? p i 

for the "heavy species" of hypernuclei. From Table Hone 
sees, however, that this relation is: i) strictly fulfilled 
only for N = Z nuclei, i.e., ^He, ^Be, fC, ^ 7 0, and fSi, 

ii) approximately correct for hypernuclei with N = Z, and 

iii) totally invalid for hypernuclei far from the stability 
line. 

Moreover, as all hypernuclei with the same A hold the 
same elementary rates (|2 . 24[) . very simple relationships 
can be established from Table U between their rates T p , 
and r„, and ratios T n / p . For instance: 



TpfiHe) 

r„(iHe) 
r„/p(iHe) 



0.8r p &Be) : 

1.6r n &Be) 

2T n/p { 9 A Be) 



0.5T P ( 9 A G), 

/9 , 



2r„ac), 

4r n/p (iC). (2.41) 



(pl\v(r)\0\) 



r 2 drji(pr)g NN v(r)g NA R 0l (V2b,r), 

(2.42) 



where g^\ and gNN are, respectively, the initial and final 
SRCs functions, one can see that the harmonic oscillator 
wave function R i(\/2&, r) is: i) picked at the origin for 
1 = 0, and ii) the integrand maxima for 1 = 1, which 
in principle should be at the distance V2& (= 2.5 fm for 
\ 2 C) from the origin, is shifted even farther because of the 
factor r 2 . This, together with the approximation (12.28)) 
for the cm. overlaps, makes the s-wave radial matrix el- 
ements large compared to the p-wave ones. Moreover, it 
could be worth mentioning that with the Wood-Saxon ra- 
dial wave functions one gets analogous results since they 
are quite similar to that of the harmonic oscillator, as 
can be seen, for instance, from 60, Figure 2-22]. 



III. COMPARISON BETWEEN EXACT AND 
sWA RESULTS 

For the hypernuclei of interest here, the approximated 
results are confronted numerically with the full calcu- 
lations in Table [TTJ This is done within the following 
framework: a) The NMWD dynamics is described by 
the 7T + K OME pot ential, with the weak coupling con- 
stants from Ref. Il8l . l24j , b) The parameter b is evalu- 
ated as in Ref. [28j, i.e.. b = l/^/hcuM-^, with Tiuj = 
45A~ 1 /3 _ 25A" 2 / 3 MeV. c) The initial and final SRCs, 
as well as the finite nucleon size effects are included in 
the same way as in our previous works [26 28, 37]. The 
results displayed in Table [TT1 clearly show that the agree- 
ment between the exact and sWA results is indeed quite 
satisfactory. In fact, the differences between them are of 
the same order of magnitude or smaller than that of the 
kinematical and nonlocality effects discussed in [27j ] . 



7 



TABLE II. Results for exact, and sWA transition rates, eval- 
uated, respectively, from Eqs. 1)2.14^ . and (|2.27p . Those for 
A Be are not shown, since they fall in between those for A He, 
and aC. 



Approx. r„ 



A ne 
A Li 



iHe 



A H 



A ^ 



r 



i 8 si 



Exact 
Exact 

sWA 
Exact 

sWA 
Exact 

sWA 
Exact 

sWA 
Exact 

sWA 
Exact 

sWA 
Exact 

sWA 
Exact 

sWA 
Exact 

sWA 
Exact 

sWA 
Exact 

sWA 
Exact 

sWA 



0.149 
0.154 
0.153 
0.265 
0.262 
0.131 
0.131 
0.208 
0.207 
0.200 
0.199 
0.205 
0.201 
0.241 
0.238 
0.340 
0.335 
0.253 
0.250 
0.279 
0.275 
0.297 
0.289 
0.311 
0.302 



0.358 
0.375 
0.369 
0.317 
0.318 
0.676 
0.637 
0.528 
0.499 
0.628 
0.594 
0.793 
0.755 
0.615 
0.582 
0.535 
0.510 
0.719 
0.689 
0.710 
0.677 
0.815 
0.760 
0.806 
0.751 



0.507 
0.529 
0.523 
0.583 
0.581 
0.807 
0.768 
0.736 
0.706 
0.828 
0.794 
0.998 
0.955 
0.855 
0.820 
0.874 
0.845 
0.972 
0.938 
0.989 
0.952 
1.112 
1.049 
1.117 
1.053 



0.417 
0.409 
0.416 
0.836 
0.824 
0.194 
0.206 
0.394 
0.414 
0.319 
0.335 
0.259 
0.266 
0.391 
0.409 
0.635 
0.657 
0.352 
0.362 
0.393 
0.406 
0.364 
0.380 
0.385 
0.401 



A more rigorous inclusion of the strong interaction in- 
gredients on the initial and final two-body states, as done 
in Refs. [U, (2|| Ho[, will modify in the same way the ex- 
act and sWA results, without affecting the conclusions 
of the present work. Namely, there is no physical rea- 
son why the mixing between states with the same total 
angular momenta and different orbital and spin angular 
momenta - induced by the SRCs, and exhibited in |4oL 
Eqs. (46) and (47)] - should influence differently the ex- 
act and sWA calculations. That this is true for the final 
state follows immediately from the fact that the sWA is 
done only on the initial state. Thus, all the discussion 
performed in Refs. [13, E3] for the final-state tensor cor- 
relation is equally valid for both calculations. The initial 
AN SRCs are less discussed in the literature. Neverthe- 
less, it was established that the phenomenological spin- 
independent correlation function [18|, Eq. (21)], which is 
the same as that used here, is a good approximation of 
the full correlation function. 

The sWA works very well for any other choice of the 
OME potential different from the model assumption of 
7r + K exchanges considered above. As one example in 



Table [TTI are also shown the results for ^ 2 C (labelled as 
A 2 C) obtained with the full tt + T) + K + p + u + K* 
OME potential. It is also evident that more realistic 
estimates of the oscillator parameter b, as the one given 
by Itonaga et al. [251 ] . would affect both calculations in 
the same manner, and would yield to the same degree of 
agreement between the exact and sWA results 0. 



160 



150 



> 

2 140 



130 



120 



0.5 

0.45 ■ 

0.4 ' 

0.35 
0.3 

0.25 i 

0.2 ' 

0.15 ' 

R, 



pi 



1.2 1.3 1.4 1.5 1.6 1.7 

b(fm) 



FIG. 1. (Color online) Single particle decay rate R p i as a 
function of the length parameter b and the liberated energy 
A P . 

One should keep in mind that in the Fermi gas model 
the A-hyperon is taken to be always in a relative s-state 
with respect to any of the nucleons within the hypernu- 
cleus. Therefore, the success of the sWA indirectly jus- 
tifies the ap plic ation of such a model to the NMWD of 
finite nuclei [39|. 

After fixing the OME potential, all Rnj depend only 
on b and Ajy. As an example, the dependence of R p \ 
on these two quantities is illustrated in Fig. [T] The 
variation of A p has a very small effect, as can be seen 
from Eq. (|2.37p . Contrary to this, the Rnj depend very 
strongly on b trough the radial wave function Roo(\/2&, r) 
in Eq. (|2.42[) . The corresponding s-shell single-particle 
decay rates Rnj are exhibited in Table Hill Finally, we 
note that by using the values listed in Tables U and IIII1 
together with Eq. (I2.27|) . we recover the sWA results 
shown in Table [TTJ 



5 The 6-values used here: 1.765, 1.781, 1.838, and 1.966 fm , for 
j^B, A 2 C, A 6 0, and A 8 Si, do not differ much from the values 
reported in Ref. |25| . which are, respectively, 1.65,1.65,1.755, 
and 1.865. 
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TABLE III. The s-shell single particle decay rates Rnj scaled 
by a factor of 10. The results for \C, and A Be are not shown 
as they are the same as those for iHe. 



A ne 
A Li 
a He 



C 



A b O 

i 7 o 



Si 



WR n0 WRm WR p0 10R 



'pi 



0.1411 
0.1310 
0.1228 
0.1162 
0.1133 
0.1107 
0.0950 
0.1038 
0.1018 
0.0860 
0.0849 



0.9470 
0.8842 
0.8334 
0.7915 
0.7732 
0.7563 
0.6557 
0.7124 
0.6997 
0.5975 
0.5905 



0.0705 
0.0655 
0.0614 
0.0581 
0.0567 
0.0553 
0.0475 
0.0519 
0.0509 
0.0430 
0.0425 



2.3620 
2.2194 
2.1026 
2.0054 
1.9626 
1.9230 
1.6845 
1.8196 
1.7894 
1.5441 
1.5272 



Tabic HTT| as well as Fig. Q] clearly show that the size 
parameter b is the most important nuclear structure pa- 
rameter for the NMWD rates r n and r p , and therefore 
the knowledge of its value for each individual hypernuclei 
could become crucial in comparing the theory with ex- 
periments. However, this not come to pass with the ratio 
r„/ p , which is mainly tailored by the OME potential. 



IV. 



CONCLUSIONS AND SUMMARY 



a function of the hypernuclear mass number, stems from 
the interplay of the increase of Fnj, and the decrease of 
Rnj, and 

(3) The ratio T n / V is almost the same for all hyper- 
nuclei that are on the stability line (N = Z), i.e., ^He, 
A Li, j^B, A 3 C, A 7 0, A 9 Si, etc. Moreover, it decreases 
when one moves toward the proton drip-line (Z > N), 
and increases when one goes toward the neutron drip- line 
(N > Z). It diminishes, for instance, by more than a fac- 
tor of 4 when going from iHe to iC, while the constituent 
Rnj rates remain the same. It might be somewhat sur- 
prising that r„ < r p even when the neutron number is 
greater than the proton number. But, as seen from Table 
IIIIl the reason for this is the dominance of 1Z p i on the 
other three single-particle decay rates. This dominance, 
in turn, comes from the dominance of the tensor ampli- 
tude d on the remaining amplitudes. The only exception 
is \H for which 1Z p i does not contribute. 

In summary, using as a tool the IPSM and the s-wave 
approximation, we have shown that the decay rates T n , 
and T p can be interrelated in a very simple way in all 
hypernuclei going from ^He up to i 9 Si. The relationships 
between them are particularly simple for the hypernuclei 
with the same mass number, as illustrated by Eq. (|2.41[) 
for the sequence iHe — >iBe— >iC. Results of this type 
are very likely valid in general, and as such they could be 
exploited to study experimentally the variations of T n , 
T p and T n / p along many similar arrays in a systematic 
way. 



The following conclusions can be drawn regardless of 
the OME potential that is used: 

(1) The sWA is sufficiently accurate, not only for qual- 
itative discussions, but also for quantitative descriptions 
of the NMWD in hypernuclei within the IPSM, when 
the SRCs are described by phenomenological correlation 
functions as done here. 

(2) The increase of transition rates T n , T p , and T nm , as 
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